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ABSTRACT 

Context. The overall properties of disks surrounding intermediate PMS stars (HAe) are not yet well constrained by current observa- 
tions. The disk inclination, which significantly affect SED modeling, is often unknown. 

Aims. We attempted to resolve the disks around CQ Tau and MWC 758, to provide accurate constraints on the disk parameters, in 
particular the temperature and surface density distribution. 

Methods. We report arcsecond resolution observations of dust and CO line emissions with the IRAM array. We also searched for the 
HCO^ J=l— >0 transition. The disk properties are derived using a standard disk model. We use the Meudon PDR code to study the 
chemistry. 

Results. The two disks share some common properties. The mean CO abundance is low despite disk temperatures above the CO 
condensation temperature. Furthermore, the CO surface density and dust opacity have different radial dependence. The CQ Tau disk 
appears warmer, and perhaps less dense than that of MWC 758. Modeling the chemistry, we find that photodissociation of CO is 
a viable mechanism to explain the low abundance. The photospheric flux is not sufficient for this: a strong UV excess is required. 
In CQ Tau, the high temperature is consistent with expectation for a PDR. The PDR model has difficulty explaining the mild tem- 
peratures obtained in MWC 758, for which a low gas-to-dust ratio is preferred. A yet unexplored alternative could be that, despite 
currently high gas temperatures, CO remains trapped in grains, as the models suggest that large grains can be cold enough to prevent 
thermal desorption of CO. The low inclination of the CQ Tau disk, ~ 30°, challenges previous interpretations given for the UX Ori - 
like luminosity variations of this star. 

Conclusions. We conclude that CO cannot be used as a simple tracer of gas-to-dust ratio, the CO abundance being affected by 
photodissociation, and grain growth. 

Key words. Stars: circumstellar matter - planetary systems: protoplanetary disks ~ individual: CQTau, MWC 758, MWC 480, 
AB Aur - Radio-lines: stars 

1. Introduction The disk around CQ Tau (HD 36910) has also been im- 

aged at different wavelengths. So far, this is one of the oldest 
It is now generally admitted that the disks surrounding inter- HAe star (-10 Myrs) surrounded by a resolved dust and gas 
, mediate mass (^ 2M0) pre-main-sequence (PMS) HAe stars disk dMannings & Sargeii^l 19971: iTesti et alJlIOOlt iDoucet etall 
(Tj • are warmer and more massive analogs of those surrounding 2006). Moreover, CQ Tau appears as a peculiar HAe star ex- 
■ ■ lower mass TTauri stars. Mannings & Sargent (1997, 2000) have hibiting an UX Ori-like variability (Natta et anil997h . MWC 
shown that disks are indeed frequent, if not ubiquitous around 758 (HD 36112) has also been barely resolved in CO by 
HAe stars. Nevertheless, only very few objects have been mod- [Mannings & Sargentj (Il997ll2000l) . Both disks appeared signifi- 
eled in detail. Accurate disk orientation, sizes, temperatures and candy weaker and smaller in CO lines than the previously stud- 
CO abundances are available fo r a few o bjects onl y: AB Aur led disks around HAe stars. 

Tp 



JPiitu et al. 2005), M WC 480 (ISimon et al. 2000; Ipi etu et all 



I2QQ7I), and HD 163296 (j lsella et al. ||2007) have been studied in ^^^^^ improve our knowledge on HAe disks, we 

CO isotopologues. Akeady one disk out of these three stands out ^ave observed CQ Tau and MWC 758 with the IRAM array 

as peculiar: the disk around AB Aur exhibits non-Keplerian mo- present here our results. The two stars were imaged in 

tions, perhaps due to its youth (Pietu et al. 2005^> |Lm et alj2006|) . 12^0 J=2^ 1 and continuum emission at 1 .3 and 3.4 mm (an- 

|HD 34382 wa^bserved at high angular resolution in ^^CO by g^i^r resolution of 1.5 and 3.4", respectively). 2.5" resolution 

iPietu etal.|(l2003D, but due to its larger distance, only limited m- images of '^CO J=1^0 in MWC 758, and upper limits on the 

formation could be obtained. HCo+ J= 1 ^0 line are also reported. 



Send offprint requests to: E.ChapiUon e-mail: Section 2 presents the Stellar properties and the observations, 

chapilloniaobs. u-bordeauxl. fr Section 3 focusses on the CO and dust modeling. We present 

* Based on observations carried out with the IRAM Plateau de Bure i" Section 4 the chemical model that we used to interpret the 

Interferometer. IRAM is supported by INSU/CNRS (France), MPG observed CO column densities. We discuss in Section 5 all the 

(Germany) and IGN (Spain). properties of the two sources before conclusion, in Section 6. 
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Table 1. Properties of HAe stars with CO disks 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Source 


Right ascension 


Declination 


Spec, type 


Effective temp.(K) 


Stellar lum.CLg) 


Distance 


CO paper 


MWC 480 


04:58:46.264 


29:50:36.86 


A4 


8460 


11.5 


140 


1,2 


AB Aur 


04:55:45.843 


30:33:04.21 


AO/Al 


10000 


52.5 


140 


3 


HD 34282 


05:16:00.491 


-09:48:35.45 


Al/AO 


9440 


29 


400 


4 


HD 163296 


17:56:21.28 


-21:57:21.9 


Al 


9300 


30 


122 


5 


CQTau 


05:35:58.485 


24:44:54.19 


A8/F2 


7200 


8-16' 


140* 


6 


MWC 758 


05:30:27.542 


25:19:57.32 


A3/A8 


8200 


11* 


140* 


6 



Columns 2&3: J2000 coordinates deduced from the fit of the 1.3 mm continuum map of the PdBI. Errors bars on the astrometry a re o f order 
< 0.07". Column s 4, 5, and 6, the spectral type, effective temperature and the stellar lumin osity are those given in 'Simon et al." l!2000) and 
iPietu et alj j2003l) for HD 34282 lvan den Ancker et all (119981) for AB Au r and MWC 758.lNatta et al. ( 2001) and Mannings & Sargent ( 1997) for 
CQ Tau. Column 8, C O interferometric papers are: 1 = |Pietu et al] j2007b . 2 = ISimon et alj j200Q) . 3 = ,PietuetaU t2005.) . 4 = ,Pietu et aU a003h , 
5 = llsella et alj i2007h . and 6 = this work. (*) For MWC 758 and CQTau, the luminosity has been scaled to the assumed distance of 140 pc. 



2. Observations 

2.1. Star properties and the distance problem 

CQ Tau and MWC 758 are two Herbig Ae stars situated near 
the edge of the Tai irus complex. Th e spectral type of M WC 758 
range s from A3e jThe et alJll994l) to A 8v (jBe skrovna ya et al 

ICreec 



1999 ): that of CQ Tau from A8 to F2 (e g. lCreech-Eakman et al 



2002h . or even F4/F5 (Mor a et al.lliOOTl) . MWC 758 has a rota- 
tion veloc ity vsin/ = 60 kms"', a nd CQ Tau v sin/ =100 km s ' 
(iBoehm & C atala 1995t iMora et al. 2001). 

Hipparcos measurements place them f ormally at 100 pc (CQ 
Tau) a nd 200 pc (MWC 758) respectively dvan den Ancker et alj 
119981) . However, the formal uncertainties on the parallaxes are 
such that the two stars could be equally placed (at the 1 .5cr level) 
at 140 pc. As the apparent position of the stars place them at the 
outskirts of the Taurus/ Auriga molecular cloud region and their 
systemic velocities are equal to the Taurus one, we assume here 
a physical association, and use D=140 pc for both sources. This 
eases the comparison with other studies of similar objects in the 
Taurus region. Furthermore, in Sec l3.4l we show that the 100 pc 
distance for CQ Tau is not compatible with our new measure- 
ments. 

Nevertheless, the effects of the distance uncertainties will 
be considered whenever appropriate in this paper Most of the 
results for the disk parameters can be simply scaled for dif- 
ferent distances us ing the appropriate scaling laws given in 
iDutrev et al.l(l2003b . 

The stellar properties corresponding to the D=140 pc dis- 
tance are given in Table [1] with those of other HAe stars having 
detailed disk studies. Note that CQ Tau being highly variable (up 
to 3 magnitudes, see Grinin et al. 2008), the values of its lumi- 
nosity quoted in the literature vary by almost an order of magni- 
tude. 



2.2. PdBI data 

The '^CO J=2-> 1 data for CQ Tau and MWC 758 were obtained 
in snapshot mode in Nov 1999 and Jan 2001, together with AB 
Aur (Pietu et al. 2005). The '-CO J=2-l data were smoothed to 
0.25 km.s"' spectral resolution for the final analysis. Baselines 
up to 170 m provided about 1.5" resolution with robust weight- 
ing. The HCO^ J=l— >0 data were obtained simultaneously: a 20 
MHz/512 channels correlator unit provided a spectral resolution 
of 0.13 km.s"'. 

The phase and secondary flux calibrators were 0415-1-379 
and 0528-1-134. The primary flux calibrator was MWC 349. 



The rms phase noise was 8° to 25° and 15° to 50° at 3.4 mm 
and 1.3 mm, respectively, which introduced position errors of 
< 0.07". The estimated seeing is about 0.2" after calibration. 
Since the observations were obtained with those of AB Aur, 
these data share a common flux scale. Thus the temperature (and 
column den sities) can be directly compared with those cited by 
iPietu et all (|2007[) for AB Aur, MWC 480 and DM Tau, as weU 
as the derived values for the dust emissivity spectral index /3. 

For MWC 758, the '2cOJ=1^0 line was observed in Dec 
1995 and Mar 1996. The longest baseline was 180 m, and pro- 
vided an angular resolution of 2.5 x 1.9" at PA 149°. 

Figure [Upresents the channel maps of '^CO J=2— >1 for CQ 
Tau and MWC 758. Figure |2] shows the '^COJ^I-^O chan- 
nel maps of MWC 758. The characteristic pattern of rotation is 
clearly apparent. No HCO^ was detected. 

We used the GILDASQ software package (lPetvl l2005) to re- 
duce the data and as a framework to implement our minimization 
technique. Table [^indicates the result of Gaussian fit to the con- 
tinuum visibilities. Elliptical Gaussian were used at 230 GHz: at 
89 and 113 GHz, the resolution is insufficient to constrain the 
disk aspect ratio and only a circular Gaussian was fitted. 

Table 2. Continuum results 



Frequency 


Flux 


Major 


Minor 


PA 


(GHz) 


(mJy) 












MWC 758 






230 


56 ± 1 


1.00 ±0.08 


0.81 ±0.10 


-10 ±20 


113 


6.7 ± 1.3 


1.1 ±0.5 






89 


3.3 ±0.4 












CQTau 






230 


162 ±2 


0.86 ± 0.04 


0.63 ± 0.04 


32 ±7 


89 


13.1 ±0.5 


0.7 ±0.3 







The error bars include only thermal noise, but not the flux scale accu- 
racy 10 %) or the effects of seeing. Here, the P.A. is that of the ellipse 
major axis. At 89 and 113 GHz, a circular Gaussian was fitted. 



3. CO and Dust Models 

The data were fitted using our standard flaring disk model, with 
power law distributions for all primary quantities (surface den- 
sity S(r) = I,o(r/roy', temperature T{r) = To(r/roy^, veloc- 
ity and scale height), following the method described in detail 



' See http://www.iram.fr/IRAMFR/GILDAS for more informa- 
tion about the GILDAS software. 
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Fig. 1. Channel maps of the CO J -2-1 hne emission and 1.3 mm continuum emission (bottom left panel) towards MWC 758 and 
CQ Tau. (a) CQ Tau Contour step is 0.075 Jy/beam or 0.56 K (2cr) for CO. Contour step is 10 mJy/beam, or 90 mK (8cr) for the 
continuum, (b) MWC 758. Contour step is 0.1 Jy/beam or 1.1 K (2cr) for CO. Contour step is 3 mJy/beam (30 mK, 2cr) for the 1.3 
mm continuum. The crosses indicate the position, orientation and aspect ratio of the disk determined from the CO analysis. The 
velocity (in km.s"') is indicated in the upper left corner. 
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Fig. 2. Channel maps of the CO J=l-0 line emission towards MWC 758. The angular resolution is 2.7 x 2.1". Contour step is 0.06 
Jy/beam or 1 .0 K (2cr) for CO. The crosses indicate the position, orientation and aspect ratio of the disk determined from the CO 
analysis. The velocity (in km.s"^) is indicated in the upper left corner. 



by iPietu et aP (l2007h . This allows the CO surface density and 
dust radial emissivity profile to be recovered completely inde- 
pendently, without assumptions on the CO abundances. An es- 
sential property of disks is that, because the CO surface den- 
sity gradient is significantly steeper than the temperature gra- 
dient ip - q > 1), it is possible, with sufficient angular res- 
olution, to distinguish (from the apparent slope of the surface 
brightness versus radial distance) the transition between an op- 
tically thick core (Ti, oc r^i) and an opti cally thin outer par t 
(Tb oc r-P or r'P-i for low J CO lines, see lDutrev et"ani2007h . 



Thus, the fitting procedure can actually reveal whether the CO 
emission is mostly optically thick or optically thin, even with a 
single transition. 

For the dust emission, we use the standard prescription for 
the mass absorption coefficient dBeck with et al. 1990..) : 



/fv(v) = 0.1(v/10'^Hz)^cm2g- 



(1) 



For the continuum data, as the angular resolution does not al- 
low to separate the optically thick core, all frequencies observed 
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with the IRAM interferometer were fitted simultaneously, using 
the temperature (fixed) derived from the CO measurements. 

The results are summarized in Table |3] The errors bars are 
derived fro m the covariance m atrix by the minimization routine. 
We refer to lPietu et al.l ( l2007h for notations (see also TableO and 
detailed description. These new observations directly prove the 
Keplerian nature of the rotation, with a best fit exponent 0.51 ± 
0.02 for both sources. 



3.1. MWC 758 

The disk outer radius is ^ 250 AU in CO and somewhat smaller 
in dust. The line and continuum emissions provide very consis- 
tent results for the disk orientation. The '^COJ=2— »1 line ap- 
pears partially optically thick, thereby allowing the derivation of 
the disk temperature in the inner regions (r < 150 AU), and of 
the CO column density in the outer regions (r > 150 AU) under 
the assumption of a single temperature power law. The param- 
eters derived from both transitions are in good agreement. The 
exponent p is not well constrained for the '^CO J=l— »0 line; we 
thus fix it to a value compatible with the '^COJ=2— >1 results. 
We also present results which both lines fitted together (col 6, 
TableO. The better sensitivity allows to derive more parameters 
such as stellar mass, temperature and exponent q. 

The difiference between the temperatures determined from 
the J=1^0 and J=2— >1 lines is in line wi th the vertical te mper- 
ature gradients found in some other stars ( Pietu et aTlbOO?!) . but 
not determined with sufficient precision to conclude. We thus 
prefer to use the temperature derived from a simultaneous fit of 
both transitions. This temperature is moderate, 30 K, with little 
dependance upon radius (exponent ^ ^ 0.3). Because the J=2— > 1 
transition is barely optically thick, the temperature is represen- 
tative of the bulk of the disk, and not only of the surface laye rs 
as in the optically thick CO disks studied bv lPietu et al.l (l2007l) . 

Because the emission is spatially resolved at 1.3 mm, the 
dust emissivity index is well constrained, with little bias due to a 
possible contribution from an optically thick core:/? = 1.0±0.15. 
If we assume that the dust and CO are at the same temperature, 
the total disk mass is 2.7 10"^ Mq (with/f(230) = 0.023 cm-g"'). 

Assuming the same distribution and excitation conditions 
than for CO, the HCO^ J=l— >0 line provides an upper limit of 
4 10'^ cm"^ for the column density of HCO+ at 100 AU, i.e. 
[i2C0]/[HC0+ ]> 10^ 

The inclination is poorly constrained (/ ^ 16°), but, ex- 
cept for the rotation velocity at 100 AU (as the data essentially 
constrain v sin /), all other parameters remain unaffected by this 
uncertainty. Besides the solution presented here, another solu- 
tion with / ^ 27° can be found. However, only the low in- 
clination solution is compatible with the stellar mass derived 
from the location of the star on a distance independent HR di- 
agram, logiLIM^) ^ 0.5 and log{T^f{) = 3.92. For / = 27°, 
the stellar mass is 0.7 ± 0.3 Mq. Reconciling this with the spec- 
tral type of the star would imply placing it at a distance of 
> 280 pc. The disk outer radius would be > 550 AU, the disk 
mass 0.08 Mq, and the temperature at 100 AU 40 K. MWC 
758 would have among the most massive disks around HAe 
stars (with HD 34282 Pietu et al. 20Q3). It would ho wever give 
log(L/M^) ^ 1.3, which from the lSiesset all ( l2000l) tracks, re- 
quire M, > 5M0. 



3.2. CQ Tau 

For CQ Tau, we obtain for the first time accurate measurements 
of the inclination i - 29 + 2° of the disk and of the orienta- 
tion of the disk axis (PA = -36 + 1°). The disk inclination is 
well constrained from both '"CO (29 + 2°) and continuum emis- 
sion (/ = 29 + 9°), although with larger error bars for the con- 
tinuum. Note that both determinations are affected by different 
instrumental effects. The dust emission is biased by the phase 
noise (seeing), while the CO-derived inclination is affected by 
bandpass calibration errors, which are negligible in our case. 
The overall agreement indicates that the seing effects remain 
small compared to the thermal noise. However, our results on the 
dis k differ from several previous determinations: for example, 
Testi et al.l (l2003h cite a major axis disk orientation PA ~ 20° 
(or -70° for the disk axis) and inclination (/ ^ 66°) very differ- 
ent from ours. Most of the previous measurements were based on 
deconvolved sizes in the image plane, either on integrated inten- 
sity maps of CO ( Mannings & Sargentll 19971) or fro m relatively 
low S /N continuum data (e.g. the 7 mm VLA data of iTesti et al.l 
1200 Ih . Our results being based on the velocity gradient and a UV 
plane analysis are much more reliable: in particular, they predict 
a stellar mass in agreement with the spectral type of CQTau. 
It is worth noting that, usin g the mid-IR camera CAMIRAS in- 
stalled on the CFHT, D oucet et al. I (l2006h have imaged the ther- 
mal emission of the dust disk at 20.5 fim and found an inclination 
angle of / = 33 ± 5° in agreement with our measurements. 

The disk of CQ Tau is somewhat smaller, 200 AU in radius, 
than that of MWC 758. Bot h continuum and CO measurements 
give the same outer radius. iDoucet et al.l (12006 ') estimated the 
disk radius to be of order of ~ 400 AU at 20.5 fim (when cor- 
rected for the distance of 140 pc). The 20 fim emission comes 
from small, warm, grains at the surface layer of a flared disk. 
On the opposite, the mm data trace large grains in the disk mid- 
plane. A difference in radial extent between these two popula- 
tions of grains is a standard feature of protoplanetary disks, as 
exemplified for example by HH 30 (iGuilloteau et al.ll2008h . 

The '^COJ=2— >1 emission appears essentially optically 
thin. Since the contribution of the optically thick core is small 
t{CO) < 1 for R > 60 AU, the disk temperature exponent q is 
poorly constrained. Best fit values are around 1, but with large 
uncertainty. We thereby adopted q = 0.5. This value is a compro- 
mise between the flat temperature distribution (q ^ 0) obtained 
in th e disk mid-plane s, as for MWC 758, or DM Tau and MWC 
480 (iPietu et a l. 2007). and the steeper laws oc curring at the disk 
surfaces, q ^ 0.65 (Guilloteau & Dutrevlll998l) . 

Even with this reduced degree of freedom, the temperature 
remains difficult to constrain. A lower limit of about 60 K can 
be obtained (in essence, because of the apparent '^COJ=2— >1 
surface brightness, taking into account that the Keplerian shear 
results in significant beam dilution), but higher values cannot 
be excluded. We adopted Tm) = 150 K; this choice has no ef- 
fect on the disk parameters, except for the derived CO and dust 
surface densities and thus disk mass and CO abundance. The 
brightness temperature being fixed, for a constant line-width and 
in the high temperature approximation, the CO surface density 
is proportional to the assumed temperature (because Th{CO) oc 
S(CO)/7'it) and the dust surface density is proportional to the in- 
verse of the temper ature (as /'^(dust) oc /c(v)Z(dust)r(dust), see 
iDartois et al] (l2()03h and their Fig.4 for details). 

If we assume that the dust and CO are at the same tem- 
perature, and the prescription for the dust emissivity in EqlT] 
which gives a-(230) - 0.035 cm-g"', the total disk mass is small, 
1.2 IO^^'Mq. Note however that this is probably a lower limit. 
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Table 3. Best parameters 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Source 


CQTau 




MWC 758 




Data 


'-CO J=2->1 


Dust 


'-CO J=2->1 


'^CO J=1^0 


'^CO 


Dust 


Vlsr (km.s"') 


6.17 ±0.04 




5.79±0.01 


5.90 ±0.02 


5.80 ± 0.02 




Orientation, PA (°) 


-36.7 ±1.3 


-36± 18 


-31 ±1 


-23 ±3 


-31 ±1 


-38 ± 7 


Inclination, i (°) 


29.3 ±1.7 


29 ±9 


18 ±6 


16 ±1 


16 ±4 


40 ±20 


Velocity law: 


^('■) = y^ooijdm 


r 










VelocityC), (km.s"') 


4.0 ±0.2 




3.6 ± 1.1 


[4.00] 


4.0 ±0.6 




Velocity exponent, v 


0.51 ±0.02 




0.51 ±0.03 


0.47 ±0.07 


0.50 ±0.02 




Stellar mass, M. ( Mq) 


1.8 ±0.2 




1.5 ± 0.7 [1.80] 


[1.80] 


1.80 ±0.5 




Z (*), (cm-0 


1.7 ±0.1 10"" 


1.7 ±0.3 10^- 


3.5 ±0.7 10"" 


1.6 ±2.4 10'" 


4.7 ±0.9 10'" 


6.0 ± 2.0 10^^ 


Sraass (*) (g-Cm-^) 




0.075 ± 0.015 








0.3 ±0.1 


Exponent p 


2.3 ±0.2 


1.3 ±0.1 


2.7 ± 0.5 


[3] 


2.9 ±0.4 


1.5 ±0.4 


Outer radius R„^t, (AU) 


200 ± 20 


200 ± 30 


300 ± 20 


230 ± 30 


270 ± 15 


180 ±40 


Temperature(*),(K) 


150 ±50 [150] 




37 ±6 


24 ±4 


30 ± 1 




Exponent q 


0.7 ± 0.5 [0.5] 




0.05 ± 0.20 


0.6 ± 0.3 


0.37 ±0.15 




SVC), (km.s-') 


0.32 ± 0.09 




0.50 ± 0.03 


0.28 ±0.10 


0.44 ± 0.02 




Scale heightf), (AU) 


22 




15 




11 




yS 


0.70 ± 0.04 


1.0 ± 0.15 



Column (1) contains the parameter name. Columns (2) and (4) indicate the parameters derived from '^CO J=2— > 1, column (5) parameters derived 
from '^CO J=l— >0, and columns (3) and (7) parameters derived from the dust emission, using the disk temperature from '^CO and the dust 
emissivity from Eq[T] Column (6) indicates the results of a sim ultaneous fi t to bo th CO lines. Note that the P.A. is that of the disk axis. SV is 
the local line width (sum of thermal -I- turbulent component see lPietu et alJ l l2007h for a description of the convention). E is the surface density 
((H-h2H2)/2) and H^^^s the mass surface density assuming a gas/dust ratio of 100. (*) values at 100 AU. Square brackets indicate fixed parameters. 
The error bars correspond to 1 cr level of uncertainties. 



because of the existenc e of vertical tempe rature gradients, as di- 
rectly demonstrated bv lPietu et all (l2006l) : if the dust is at 50 K, 
the disk mass increases to 3.6 10"^]V[o.We stress that our choice 
of temperature maximizes the apparent CO abundance, as this 
scales as T^. 

Our result on y6 = 0.70 ± 0.04 agrees very well with the 
j6 - 0.6 + 0.1 determined bv lTesti et al. (2003) using a broader 
frequency range including VLA observations at 7 mm. 

Assuming the same distribution and excitation conditions 
than for CO, the HCO^ J=l— »0 line provides an upper limit of 
4 10'^ cm"^ for the column density of HCO+ at 100 AU, i.e. 
[12CO]/[HCO+]>4101 



3.3. Prediction for tiie CO J=3-2 Line 

The '^CO J: =3^2 in CQ Tau and MWC 758 was detected by 
iDent et all ( l2005h using the JCIVIT. Our best fit model can 
be used to predict the intensity of this line. This is given in 
Table m also for 1VIWC480 using the model parameters from 
[Pietu et all (|2007|) . The agreement is excellent for MWC 758 and 
MWC 480, but our prediction is a factor 2 too high for CQ Tau. 
The discrepancy can be somewhat reduced (by about 30 %) by 
using a temperature of 70 K for CQ Tau. The remaining dif- 
ference may be due to noise, as the CQ Tau CO J=3-2 line is 
detected with a S/N of only 6 bv Dent et al.l(l2005b . The spectral 
baseline uncertainty on the JCMT data can also play a signifi- 
cant role at this sensitivity level. Nevertheless, this result sug- 
gests that the temperature in CQ Tau may not be as high as we 
estimated. It is also possible that the '^CO J=3— >2 line is some- 
what sub-thermally excited, as its higher opacity makes it more 
sensitive to the upper layers than the J=2— >1 transition. 



Table 4. '^CO J=3^2 comparison 

Source ''^C0J=3— »2 integrated line flux (Jy kms"') 
name observed predicted 
CQTau 6±1 14(#)-11(*) 
MWC758 15 ± 2 15 

MWC480 52 ±1.1 47 

'^CO J=3->2 observation from lDent et al.l j2005 '). The conversion fac- 
tor from the main beam brightness temperature Tn,b(K) to flux density is 
19 Jy/K for the 22" beam at 345 GHz. (#) for T=150 K,(*) for T=70 K. 



3.4. Distance effects 

As mentioned in Sec|2Tl CQ Tau could be at 100 pc and MWC 
758 at 200 pc, instead of our assumed D=140 pc. The surface 
density scales as D'', the temperature as D', the disk mass as D^, 
and the st ellar mass, when d erived from the Keplerian rotation, 
as D (see Dutrev et al.ir2003!) . With the values of p and q from 
Table[3] the impact of the distance ambiguity remains limited to 
a factor 2 on the disk masses. However, the difference in disk 
mass between CQ Tau and MWC 758 becomes affected by a 
factor 4. 

The CO abundance (at 100 AU) will scale as j:>''^co)-p('i"''K It 
is thus little affected by the distance uncertainty (less than 50%). 

Figure |3] shows the location of CQ Tau and MWC 758 on 
a distance independent HR diagram, as used by 'Sim on et al.l 
(2000), using the evolutionary tracks of Siess et al. (200(1. The 
location of CQ Tau in this diagram is a further argument for the 
distance of 140 pc for this object: at 100 pc, the kinematic data 
would require a mass of 1.3 M©, which is clearly incompati- 
ble with the effective temperature and L/M^, which indicate a 
~ 1.8 Mq star. We point out that, with these new results, CQTau 
could actually be younger than initially thought (5 Myr), and 
even younger than MWC 758. 
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Fig. 3. Distance-independent HR di agram (L/M^ ver sus Teff). 
Theoretical evolutionary tracks from Siess et alj (l2000h and Z - 
0.02 metallicity for stellar mass between 1 and 4 M© (from left 
to right: 1 to 2 by 0.1 Mq, and then 2.2, 2.5, 2.7, 3.0, 3.5, and 
4.0 Mq). Knots on the curve begin at and are spaced by 1 Myr 
CQ Tau is the leftmost cross, MWC 758 is the rightmost one. 

3.5. Comparison between sources 

Both sources have many characteristics in common: 

1 . Partially optically thin CO line emission, in contrast to most 
previously studied disks. 

2. Different surface density profiles for CO and dust, with the 
CO surface density falling faster (p ^ 2.5 - 3) than the 
dust opacity (p ^ 1.5). This property is apparently shared 
with all other observed sources so far, with the exception of 
HD 163296 which displays a relati vely shallow CO s urface 
density profile (p ^ I, according to llsella et alj (|2007|) ). 

3. Similar values for the dust emissivity exponent = 0.7-1 . In 
both sources, we have resolved the continuum emission, so 
this value is free of the bias due to a possible optically thick 
core. For CQ Tau. lTesti et al.l ( l2003h also found that the free- 
free contribution is negligible. These values are smal ler than 
those found for MWC 480 and ABAur, ^ 1.4 (PietuetalJ 
l2006Ll2007h . 

4. Similar lower limit for the [i2cO]/[HCO+] ratio, ^ 4000, 
corresponding to ['^CO]/[HCO^] > 100 when taking into 
account a reasonable (average) CO isotopologue ratio. This 
limit is co nsistent wi t h wha t has been measured in other 
sources bv lPietu et all (l2007h . 

5. An essential result, which remains whatever the distance, is 
that the apparent CO abundance is very low in both sources, 
less than 10 corresponding to a CO depletion > 10 com- 
pared to typical molecular clouds jOhishi et alJll992l) . 

6. Both sources have also relatively small disk masses, less than 
or about 5 10"^ Mq. Because of the distance problem, and 
dust temperature uncertainty for CQ Tau, disk masses remain 
the least accurate parameters. In the most extreme case, the 
MWC 758 disk could be up to 8 times more massive than the 
CQ Tau one. 

Although the two sources are very similar, CQ Tau stands 
out among the sources studied so far because it accumulates a 
number of peculiarities. It is at least as warm as AB Aur, the 
warmest source detected in CO J=2-l so far, and only the sec- 
ond source in which '^CO J=2— >1 is optically thin, after BP Tau 



dDutrev et alj|2003h . Like BP Tau, the CO and dust outer radii 
are identical within the current error bars. This is an important 
difference with most other sources, either TTauri stars or HAe 
stars, for which the outer radius of the bulk of the dust emission 
is in general smaller than the CO outer radius (see the T Tauri 
stars DM Tau, LkCal5 and the HAe sta rs ABAur, MWC 480 
(iPietu etalJl2006ll2007h . and HD 163296 (llsella et al.ll2607h '). 

3.6. CO abundance: an anomalous gas to dust ratio ? 

For both sources, one of the surprising findings is the very low 
apparent CO abundance, ^ lO"*" (within a factor 3, as this value 
decreases with radius), if a standard gas to dust ratio of 100 is 
assumed, i.e. an apparent depletion of ^ 100. Note that this re- 
sult is completely distance independent. Since the temperature 
of these two disks is large, depletion of CO due to sticking on 
grains cannot be invoked in a simple way. The situa tion here is 
simila r to that found for the T Tauri star BP Tau bv iDutrev et al.l 
(I2003h . 

This suggests either i) that the gas (H and H2) to dust ra- 
tio is actually much lower than our assumed value of 100, be- 
cause we are observing disks in the process of dissipating their 
gaseous content, or that the dust distribution has evolved (grain 
growth and settlement along the mid plane) either ii) by simple 
coagulation and size redistribution, and/or Hi) by mantle accre- 
tion, taking off elements from the gas phase. In all cases, the 
CO/dust ratio is no longer a direct reflection of the gas to dust 
ratio, because of changes in the H2 shielding of the CO photo- 
dis sociation lines. 

iJonkheid et alj (l2007l) have studied the chemistry of HAe 
disks, and find that for very low masses of small grains, ^ 
10"^ Mq, the CO/H2 ratio can be of order of a few 10"^ (see their 
models B4/BL4). Thus, our low observed CO column densities 
might be the result of grain growth and/or efficient dust settling. 
However, only in the lowest mass disks are the predicted col- 
umn densities as low as we observe, a few 10'^ (their Fig. 8). 
Asljonkheid et al. (2007) do not present calculations with grain 
growth for disk masses in the range of what we observe, we have 
performed new chemical models which are described in the next 
section. 



4. The Chemical Models 



We us e the PDR code from th e Meudon Group (" Le Bourlot et al.l 
'19 931: iLePetitetalJ I2OO6I) . with the modifications from 
Hil v-Blant et"al] ( |2008|) for the grain size distribution. The model 
is a one-dimensional stationary plane-parallel slab of gas and 
dust illuminated by an ultraviolet (UV) radiation field. The ra- 
diative transfer in the UV (which takes into account the self- 
shielding of H, H2 and CO lines and absorption in the con- 
tinuum by dust grains), the molecular abundances and option- 
ally the thermal balance are calculated iteratively at each point 
in the "cloud". The ch emical network is similar to that of 
iGoicoechea et al.l (l2006l) . No freeze-out onto grains is consid- 
ered: this assumption will be discussed in Sec |5.4| 

In order to study the chemical and physical effects of grain 
growth on disks, in particular the UV p enetration, the PDR code 



the UV pe 

was modified bv lHilv-Blant et afl (l2008i) to introduce a standard 
power law grain size distribution n(a) oc a"'' with a+ and a_ be- 
ing the maximum and minimum cutoff radii, respectively. This 
modification affects the UV extinction curve, the chemistry (H2 
formation) and the thermal balance and is self-consistently in- 
troduced. The resulting extinction curve is calculated using the 
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Table 5. Modeling parameter values 



Parameters 


values in calculations 


gid 


1, 10, and 100 


a+ 


1 mm, 10, 1, and 0.1 yum 


X (at 100 AU) 


10^10^ 5 10\ 10"*, plus a modified shape 


AV(km.s-') 


0.5, 1.0, and 2.0 



gId is the Gas to Dust ratio, a+ the maximum grain radius, x the UV 
scale factor over the Draine interstellar field, and A V the local line width 
for H2. See Section 4.1 for details. 

Mie theory for homogeneous isotropic spherical particles and 
presented in AppendixlAl 

To have a two-dimensional molecular distribution we com- 
pute the model at different radii. The output of this 1 + lD model 
is the vertical distribution of molecular abundance calculated at 
different radii. 

4. 1 . Parameters of the Chemical Model 

As an input we can impose physical conditions relevant from the 
disk structure (i.e. temperature and vertical density laws as de- 
rived from TableO at each radius. Since the disk masses are the 
most poorly constrained parameters, the models are not intended 
to fit any particular object, but are to be used to understand the 
major effects of the model parameters on the CO abundance, 
and more specifically, to seek what conditions and physical pro- 
cesses prevent CO to reach high surface densities. Two different 
disk density profiles were used. The "low mass" case correspond 
to a dust mass (in the inner 200 AU) of 11 10^'' Mq, the "high 
mass" case of 33 10"^ M©, both with p = 1.3 and a scale height 
of 22 AU at 100 AU. 

We investigate several gas to dust ratio (g/d) under various 
UV field conditions, and different maximum grain size a+. A 
wide range of parameters has been explored (see Table |5]). The 
whole ensemble of results is presented in Appendix lAl IB] and ICl 

UV field To estimate the UV field, we assumed UV is scat- 
tered by the dust: half towards the mid-plane and half towards 
the exterior. We have checked that the opacity toward the star is 
strong enough to justify this hypothesis, as the amount of small 
grains remains sufficient even when increasing a+. The shape of 
th e UV field is assumed to foll ow a Draine field (as provided 
bv lSternberg & Dalgarnol ( ll995l) ). With the above scattering as- 
sumption, it is described by a scale factor (0) at 100 AU, and 
decreases as l/R^. As the actual UV fields from the stars are 
unknown, several simulations have been done with different x 
factors (lO"* and lO'*). As expected, a decrease of the incident 
UV field implies a decrease of the photodissociation rates, and 
makes CO more abundant (Figure |BTT] and |R4] in Appendix B). 

Thermal Balance Models can be calculated with and with- 
out thermal balance . The dust temperature is derived following 
iBurton et al.l (Il990h . With thermal balance, this dust tempera- 
ture is used to compute the energy transfer between the gas and 
the grains as described in Burke & Hollenbach ( 1983). The gas 
temperature is then cal culated, taking into account heating and 
cooling processes (see iLe Petit et all (12006,) ). Without thermal 
balance, the gas temperature is taken from the CO observations. 
The thermal balance has little effects on the CO abundances, ex- 
cept at very low densities. Provided that the assumed gas and 
dust temperatures are close enough to those derived in the disk 

^ is the ratio of the local UV field over the Draine interstellar UV 
field 



mid-plane from the thermal balance, the results on the disk struc- 
ture and CO column densities are rather similar This is clearly 
seen by comparing Fig |6] and Fig lA.3l 

Grain Distribution The composition of the dust is 50% of 
silicates and 50% of graphite with y - -3.5 and «_ =3nm in 
all calculations. The minimum grain radius is small enough 
to properly represent the photoelectric process, the UV extinc- 
tion curve and the formation of H2 on grain surface. Keeping the 
dust mass constant, we have simulated grain growth by varying 
the maximum radius «+. The amount of small grains is there- 
fore reduced to the benefit of large ones and the extinction curve 
is clearly modified. With y = -3.5, the UV opacity scales as 
1/ -^fl7forfl+ > lOyum. 

We present the various models in Appendix A, with 
fl+varying from 0. 1 fim to 1 mm. The main result is a significant 
decrease of the opacity towards the mid-plane and therefore an 
extension of the photo-dissociation layer when a+increases. As 
a consequence, the transition between the region where Carbon 
is predominantly atomic or ionized, and the region where CO 
dominates appear closer to the disk mid-plane (figure IA.2I ) and 
the CO column density is reduced ( |A.3t . 

g/d ratio We modify the ratio by decreasing the gas mass 
(H) while keeping the dust mass constant to simulate the dissi- 
pation of the gas. We performed calculations for g/d = 100 and 
g/d - 10 with several grain size distributions. The main effect of 
decreasing the gas to dust ratio is a modification of the spectral 
distribution of the radiative energy density in the disk (see figure 
nil. The gas being less abundant, the absorption in the lines of CO 
and H2 is less efficient and the transition H - H2 occurs closer to 
the mid-plane, and the disk is more ionized (e.g. for a+ - lyum, 
C^ is the dominant form of carbon in the mid-plane for radius 
greater than 200 AU). 

Turbulence We also check the effect of the intrinsic line 
width (turbulence) on the photo-dissociation line shielding. We 
have made three runs with three different values of the Doppler 
width of the H2 lines (0.5 - 1.0 - 2.0 kms ') for a+ - lyum and 
1 mm with;^^ = 10"*. The stronger effects (obtained for a+ - \ 
mm) are presented in Appendix |C] but the effect on CO is neg- 
ligible, essentially because the CO pre-dissociation intrinsic line 
widths are much larger. 

4.2. Chemical Models and comparison with the CQ Tau Disk 

We only discuss in this section the results relevant for the CQ Tau 
case. The MWC 758 case is treated in the discussion by compar- 
ison to this example. More models are given in Appendix lAl iBl 
and|C] We assume the physical and chemical conditions derived 
from the CO analysis (see TableO. Our goal is to reproduce the 
low CO column density which is observed, and corresponds (for 
a normal gas-to-dust ratio) to a CO depletion factor of ~ 100 or 
an abundance CO/H2 ^10"^ 

Figures|5]and|6|shows the abundances with respect to H-1-2H2 
and the column densities obtained at radii 100, 200 and 300 AU 
under a UV field of x - 10^. We assume two different values 
for a+{lfim and 1 mm) and for g/d (10 or 100) and the ther- 
mal balance is calculated. The gas temperature goes up to sev- 
eral thousand Kelvin in the disk atmosphere. In the mid-plane 
in the small grain case, the gas temperature is of the same or- 
der of magnitude as the dust temperature because the density is 
high enough to have thermal coupling by collisions. The shape 
of the CO abundance distribution is clearly affected by reactions 
at high temperature with OH leading to a secondary peak not 
located in the disk mid-plane. Although this peak is quite high 
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Fig. 5. Vertical distribution through the disk of the abundance of H, H2, C^, C and CO and gas (dashed line) and dust (dotted 
line) temperature at the radii 100, 200 and 300 AU for the models with standard UV field = 10"* at 100 AU), a+ - I /jm (top), 
a+ = 1 mm (bottom) and with the thermal balance calculated. Dust temperatures are plotted for the extreme grain sizes (a+and a-) 
and an intermediate value C'^^"' ). In the big grains case (bottom), the dust disk is so optically thin that we never reach Ay = 1. Note 
that the vertical and radial scales are the same in all figures. 



in abundance, it occurs when the density is low and the column 
density is very weakly afifected (Figure|6]l. 
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Fig. 6. Radial distribution of the column density of C^, C and CO for the model with standard UV field, the thermal balance 
calculated, a+ = 1 /im (left) and a+ = I mm (right). 
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Fig. 4. Energy density between 900 and 1200 A for gas/dust - 
10 (top) and 100 at R=200 AU and Av = 0.01 with a+ = Ifim, 
10^ AV=lkms-i 



4.3. Chemical uncertainties and comparison witin otiner 
results 

iJonkheid et al.l (l2007h have also calculated the chemistry and gas 
temperature of evolving Herbig Ae disks. Their model differs 
in several ways from the above approach. Contrary to us (see 
above), they have modified the gas to dust ratio by decreasing 
the mass of the small grains (i.e. interstellar grains) in order to 
mimic the dust growth and/or settling. The UV transfer is per- 



formed with a 2-D code dvan Zadelhoff et aLll2003h whereas in 
our case this is a 1-D approach. But the self-shielding of H2 and 
CO are calculated assuming a constant abundance whereas we 
compute it explicitly at e ach point of the "cloud ". T heir chemical 
network is described in iJonkheid et alj (|2004|) and iJansen et al.l 
(11995 *). It also incorporates chemical reactions with PAHs which 
are not considered in our model. 

Indeed, we also find that CO becomes underabundant if and 
only if photodissociation is important. This requires a reduc- 
tion in the UV opacity which is obtained by grain growth to 
1 mm, or larger Note that dust sedimentation towards the disk 
mid-plane should not play a significant role, as the small dust 
grains which are responsible for the UV opacity remain coupled 
to the gas. However, contrary to i Jonkheid et alj (1200 7), we find 
that Carbon is predominantly in its ionized form C^ when CO is 
photo-dissociated. 



5. Discussion 

5.1. CO Abundance: Photodissociation with grain growth 

Figure |6] (right) suggests that a case with a+> 1 mm and g/d - 
100 can explain CO column densities of order 10'^ cm"^ around 
200 AU. Figure |6] (left) allows us to conclude that the case with 
a+ - I jjm cannot explain the observed CO column densities, 
even with a low g/d ratio. For CQ Tau this is in agreement with 
the spectral index we measure (J3 - 0.7) which in dicates that 
significant grain growth has occurred as also found bv lTesti et al.l 
(.200 3). Grain growth is also supported by the (3-1 value found 
for MWC 758. 

This result suggests that grain growth, or more precisely the 
enhancement of the UV penetration resulting from grain growth 
is the dominant process explaining the measured CO column 
density. Note that we have assumed that the dust size distribu- 
tion follows dn{d) cx with y = -3.5. Th e conclusion w ill not 
significantly change for realistic slopes y dPraind (12006) indi- 
cates -y - -3.2 for /3 - 0.7). As the extinction curve is governed 
by the amount of small grains, it is reasonable to conclude that 
any solution which guarantees the same amount of small grains, 
leading to a similar extinction curve, would be sufficient to ex- 
plain the CO column density without changing g/ii. Note that for 
-y = -3.5, the mass within grains smaller than a is about yfafal 
times the total dust mass, i.e. 3 % of the dust is in grains smaller 
than Ijum for a+ - I mm. 
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It is worth pointing that the photodissociation offers a natural 
explanation for the steeper surface density law of CO compared 
to dust. As UV penetration is more efficient in the outer parts of 
the disk, CO is more heavily dissociated, and the CO abundance 
decreases with radius (see Fig|6]l, resulting in a larger slope p for 
CO than for dust and Hydrogen, in qualitative agreement with 
the observations. However, the p values predicted by the models 
(4-6) seem larger than observed (2.5-3). 

The [^^CO]/[HCO^] ratio can also be used as a secondary 
diagnostic of the physical conditions. For large grains, a low gjd 
ratio predicts a value around 1000 for this ratio, and much higher 
values for gjd - 100. Although this is relatively inconclusive 
compared to our current limit (4000), longer integration times 
on HCO^ could bring a useful diagnostic. 

Note however that the CO column density is very sensitive 
to the UV penetration (and thus to the disk mass). We computed 
models twice less massive than those presented in Figs|5]|6] In 
such cases, C^ is the dominant form of Carbon throughout the 
disk, and the CO column densities become much smaller than 
observed. Similar results can be obtained using a larger maxi- 
mum grain size. 
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Fig. 7. Flux at 100 AU of FUSE and lUE observations MWC 758 
(points), an A3V star according from the Kurucz atlas, a black 
body with the same temperature as an A3 star and several scaled 
Draine fields. The 10^ one is calculated to have the same inte- 
grated intensity between 912 and 2400 A as the observations. 
The star is assumed to be at 140 pc. 



5.2. The UV problem 

A significant uncertainty is linked to the knowledge of the UV 
field and to its vertical diffusion in the disk. Comparing Figure 
IB. 3 1 with Figure |A3] in the a+- 1mm case shows that decreas- 
ing the UV field to;^^ = 10^^ essentially brings back the Carbon 
into CO and increases the CO column densities by a factor 100 
in both cases (g/d = 10 and 100). It is important to stress 
that what matters here is the number of available photons be- 
tween 900 and 1200 A where all the important photodissocia- 
tion processes occur. Such large values of the UV field cannot 
be provided by the relatively cool stellar photospheres, and must 
come from the UV excess due to accretion, see Appendix for 
details. On the other hand, the heating will be significantly af- 
fected by the mean intensity up to 2000 A, which is dominated 
by the stellar photosphere. In this respect, note that many au- 
thors specify t he strength of the radiation field in reference to 
iHabingl (Il968l) . i.e. by the integral of the UV flux between 912 
and 2400 A (=5.166eV) or between 912 and 2050 A (=6e V) 
(I van Dishoeck et aT]|2006l) . As the shape of the UV spectra from 
T Tauri or Herbig Ae stars differs quite significantly from the 
interstellar UV field between 900 and 2400 A, this can lead to 
some inappropriate results. 

For example, the FUSE and lUE spectra for MWC 758 at 100 
AU from the star (Figure |7]) are relatively well fit, between 900 
and 1200 A, by a Draine field scaled up by;\f ^ lO'* (the apparent 
best fit is for a few 10^, but cor rection for the star extinction, 
E(B-V)=0.07. lMartin-Za"aill2005t must be added to the curve 
presented in Fig|7]l. On the other end, a scale factor of x - 10^ 
would be needed to match the integrated intensity between 912 
and 2400 A. This scale factor would overestimate by a factor 10 
the available UV flux in the 912 - 1200 A domain relevant for 
photodissociation. 

For CQ Tau, the situation is less clear because of its variabil- 
ity. From the lUE spectrum, the CQ Tau UV flux is 10-30 times 
lower than that of MWC 758 (see lGradv e"tani2005l their Figure 
5). However, the star may have been observed during a minimum 
by lUE. We performed another calculation with big grains and a 
lower UV field (x = 10^). The results for g/d - 10 are in good 
agreement with the observations: the CO column density is about 
lO'^cm"^ at 200 AU (FiglSjl and the gas temperature determined 
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Fig. 8. Radial distribution of the column density of C^, C and 
CO for the model with a+= 1mm andx - 10^. 



from the thermal balance is about 60 K. The ['2cO]/[HCO+] ra- 
tio is consistent with our limit. Thus a solution with gas dispersal 
cannot be excluded if the actual UV field is indeed relatively low 
in CQ Tau. 

5.3. MWC 758 and the Temperature problem 

The observed CO temperature and that calculated from the (ap- 
proximate) thermal balance are in reasonable agreement for CQ 
Tau, around 60 - 100 K. MWC 758 poses a different challenge, 
as the observed gas temperature is much lower, about 30 K, 
while the estimated CO column density is similar. 

Reducing the UV flux inside the disk is the simplest way 
to obtain a lower temperature. For a similar incident UV flux, 
increasing the attenuation through the disk, i.e. having smaller 
grains, will lower the temperatures. However, this will also re- 
duce photodissociation and produce more CO. Indeed, Fig|5] 
shows that low gas temperatures (in the mid-plane) can be ob- 
tained only with small (maximum size Ipm) grains, in which 
case the CO abundance is high. For all cases compatible with 
low CO column densities, the gas temperature is high (> 50 K). 
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In fact, in the PDR interpretation, to have a similar CO surface 
density, the UV fields in the disk must be similar in both sources. 

What can explain the temperature difference between the two 
objects? The thermal balance of the gas is controlled by several 
mechanisms. In the upper layers, the energy release due to H2 
formation dominates the heating, and the gas cools by radia- 
tion in the fine structure lines of C and O. Deeper in the disk, 
photoelectric effect on dust grains become important (and some- 
times dominant), and then the energy deposition from cosmic 
rays (which is largely transported via the chemical reactions). 
CO line cooling becomes important. At still higher densities, the 
gas-grain coupling become dominant: this can be either a heating 
or a cooling mechanism, depending on whether dust is warmer 
than the gas or vice versa. 

The similarities and differences between CQ Tau and 
MWC 758 are 

- CQ Tau apparently has a lower UV field than MWC 758, by 
a factor 10 to 30 (see Sec|5^ 

- The UV attenuation is most likely larger in MWC 758 than 
in CQ Tau, because of the larger density and the (slightly) 
smaller grains (as its fi value is larger). 

- The smaller grains in MWC 758 lead to a higher dust tem- 
perature, although the effect is expected to be small (e.g. 
iKruegel & Sieb enmorgen' 1994V 

- The larger densities in MWC 758 result in a more efficient 
gas-grain coupling, as the coupling constant for thermal ex- 
change 

«dust«gas gid 
A oc oc (2) 

-\/a-a+ -\/a-a+ 

(see lBurke & Hollenbachll983h . where we assume the grain 
size distribution has an exponent y = -3.5. 

- The above equation also indicates that smaller grains in 
MWC 758 also result in a more efficient gas-grain coupling. 

Can these differences explain the lower gas temperature in 
MWC 758? In addition to the above effects, dust sedimentation 
may play a role. If dust grains are more sedimented in CQ Tau, 
the gas-grain coupling will be less efficient, and this may pre- 
serve a higher temperature of the gas above the disk mid-plane. 

The magnitude of these effects is difficult to evaluate with- 
out detailed modeling. Figure |9] indicates the heating/cooling 
rates for a disk with a (dust) surface density compatible with 
that of MWC 758 as computed in our models. In this figure, the 
"chemical" heating rate is the total energy released by chemical 
reactions: one of the prime energy source for these is actually 
cosmic-rays through dissociative recombination ion-neutral re- 
actions, the primary ions being pro duced by cosrnic-ray s. Note 
that it is probably an upper-limit (see lLe Petit et alj2006l section 
6.1.5). The direct heating by cosmic-rays is too low to appear 
here. 

Taken at face value, the results of the thermal balance study 
would favor the low g/d solution (i.e. g/d - 10), which results 
in a somewhat more efficient cooling. Reducing even further 
the g/d ratio (g/d = 1) results in higher gas temperature be- 
cause the gas is then mostly heated by H2 formation and pho- 
toelectric effect on grains. However, there are several signifi- 
cant uncertainties in our modeling procedure. For example, the 
amount of small grains which controls the photo-electric heating 
dBakes & Tielensl[l994 ') is poorly constrained, so the efficiency 
of this process may be overestimated in our model. Similarly, 
the gas-grain coupling is dominated by the small grains, because 
of their larger cross section and higher temperature, and remains 



also rather uncertain. Furthermore, we use a very crude approx- 
imation for the diffusion of UV photons towards the disk plane, 
and even the unattenuated UV flux is uncertain (see Sec l5.2l i. 
Note that the photodissociation is totally dominated by the UV 
excess, while the stellar UV flux will play a role in the heating 
processes. 

In fact, the situation of CQ Tau and MWC 758 is comparable 
to that of BP Tau, where Dutrev et al. (2003) have found warm, 
optically thin CO, with a very low "apparent" abundance (of or- 
der 10"^). It is thus tempting to interpret BP Tau in the same 
framework of CO photodissociation by the UV. This possibil- 
ity is further supported by the fact that BP Tau is a very active 
T Tauri star, with a large accretion rate (2 lO^^Mo/yr, see e.g. 
iBertout et alJlT988l:lGullbring et alJlT996l) . 

5.4. An alternative: Grain growtii by ice mantle accretion 

We have not explored so far the possibility that Carbon (and 
Oxygen too) nuclei are locked into grains. The gas temperature 
determines whether gaseous components can stick onto grains, 
while the dust temperature determines whether ices can evapo- 
rate. In both sources, the gas is warm enough to avoid new CO 
molecules to stick on grains. However, a fraction of C and O 
nuclei can nevertheless remain trapped in grains, as discussed 
below. 

This can happen if Carbon has been removed from the gas 
phase through the building of ice mantles onto grains. Such a 
process does not change the dust mass significantly (at most a 
factor 2 with standard elemental abundances), and changes even 
less the grain sizes. However, it removes selectively from the gas 
heavy nuclei with respect to Hydrogen. It is then conceivable to 
have a low CO / dust ratio. The difficulty of such a possibility is 
the thermal history of the dust and gas during the disk evolution. 
The dust must have been cold enough to accrete CO, which is the 
main reservoir of Carbon, and not released later, either thermally 
or by photo-desorption due to the enhanced UV field. 

This can happen if the CO is converted onto grains into other 
species (such as H2CO or CH3OH) with larger desorption tem- 
perature. 

In this respect, it is worth noting that in our simulations, the 
temperatures of the big grains are low enough (around 10 K) to 
prevent CO from being thermally evaporated. In a dynamical 
process where the small grains build up ice mantles, and then 
grow by coagulation to form larger and larger grains, the bulk of 
the CO trapped on grains will end up in the large grains, which 
contains most of the mass. Hence, despite the apparent "large" 
gas temperature traced by CO, a substantial fraction of C and 
O may be actually depleted onto grains. In this explanation, 
the apparently low CO to dust ratio would not imply a low g/d 
ratio, but strong C and O depletion like in disks around T Tau 
stars. The achieved depletion depend on the past thermal history 
but also critically on the grain growth processes. The maximal 
depletion depends on the dust mass fraction in grains which are 
large enough to remain below the CO desorption temperature. 

We have computed some models with such depleted elemen- 
tal abundances. As expected, the CO abundance reflects the ele- 
mental depletion (see FigurefTOl). Note that the gas remains warm 
even in the "dense" case (factor 3 in density compared to the pre- 
vious models). 

A possible way to check this hypothesis would be a direct 
measurement of the mm dust temperat ure. The onl y such avail- 
able measurement is for MWC 480 by iPietu et alj (j2006), who 
indeed found a very low dust temperature in this source, much 
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Fig. 9. Heating (left) and cooling (right) process for the models with small grains (top) and big grains (bottom) and a gid ratio of 
100 in erg cm'-' s'' . H2-gr is the heating due to H2 formation on grains, ph-gr is the photoelectric effect on grain, photon is the direct 
heating of the gas by the photons and chem is the heating due to chemistry reactions (see text for details). Squares: gas temperature, 
triangles: upper and lower limit of dust temperature. 



lower than found from the SED fitting, and perhaps as low as 10 
K near 100 AU. Although this can be attributed to vertical tem- 
perature gradients, the temperature dependence on the grain size 
may also play a role in explaining such a large difference. 

5.5. The pending UX Ori Problem for CQ Tau 

The large photometric variations observed in CQ Tau classify it 
as an UXOr variable, and have been interpreted as due to vari- 
able extinction. The accurate determination of the inclination of 
CQ Tau challenges this interpretation since UX Ori phenomenon 
can only be observed if the disk incUnatio n is larger than ~ 45" 
(iNatta et al.l[T997l:lNatta & Whitnevll2000l ). 

lEisner et"aLl ( |2004|) have measured the distribution of the 
NIR emission from the disks of CQ Tau and MWC 758 using 
the PTI interferometer. For CQ Tau, they find a position angle 
for the disk axisE of 14 ± 6° and / = 48 ± 5° (for MWC 758, they 
find PA= 40 ± 6° and ; = 33 + 4°). Taken at face value, this in- 
dicates that the orientation and inclination of the two inner disks 
differ markedly from that of the outer disks, sampled by CO, the 
axes being separated by 30 ± 7°. One possibility would be that 
the inner disk is warped due to dynamical interactions with (yet 
undetected) i nner bodies, as in the case of the more evolved disk 
of jS Pictoris (Mouillet et al.|[T997b . Such large warps have never 



^ lEisner et give the PA of the emission major axis whereas 

we indicate the PA of the disk axis 



been observed so far. In the case of t he spectroscopic binary UZ 
Tau-E, for which lsTmon et aL I (120001) suggested a mis-alignment 
between the orbit and the surrounding disk by a similar amount, 
a recent re-evaluation of the binary orbit (Prato et al. 2002) con- 
cludes that the orbit and disk planes are aligned to within 4 de- 
grees. In all confirmed cases, disk warps are small: only a few 
degrees in /3 Pic for example. 

Moreover, the inner disk results are based on very few inter- 
ferometric data points: a confirmation of theses values would be 
required before one can seriously invoke such large disk warps. 

6. Summary 

Using the IRAM array, we have observed in continuum and CO 
J=2-l and J=l-0 the disk surrounding MWC758 and in CO J=2- 
1 that of CQ Tau. We also got upper limits in HCO^ J=l-0 in 
both sources. We detected the dust and gas emission and were 
able to derive their physical conditions applying a standard disk 
model. In both cases, we found two relatively small CO and dust 
disks compared to those observed and analyzed so far around 
other Herbig Ae stars such as AB Aur, MWC 480 or HD 34282. 
We also used the PDR model developed by the Meudon Group 
in order to qualitatively interpret the observed CO abundances. 
Our main results are: 

- The two disks have partly optically thin CO J=2-l line and 
high CO depletion factor with respect to H2 of the order 
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Fig. 10. Vertical distribution through the disk of the abundance of H, H2, C^, C and CO and gas (dashed line) and dust (dotted line) 
temperature at the radii 100, 200 and 300 AU for the models with standard UV field {x=lO*), a+- 1mm, g/d-lQQ in the dense 
case (MWC758 like). Top: all the elemental abundances (except H and He) are depleted by a factor 10. Bottom: standard elemental 
abundances 



~ 100 - 200. We also found a similar lower limit on the 
ratio [12C0]/[HC0+] ^ 4000. In both cases, their CO sur- 
face densities fall down faster than the dust surface density. 
The dust emissivity spectral index is of the same order with 
jS ^ 0.7 - 1, implying grain sizes up to a few mm. 

- The main difference between the two disks is the tempera- 
ture. The disk orbiting CQ Tau is the hottest one observed so 
far with Tioo ^ 150 + 50 K, while the MWC 758 disk is not 
as warm, with Tioo - 30 + 1 K. 

- Although the disk masses are uncertain, due to distance am- 
biguities and, in CQ Tau, to limited knowledge of the temper- 
ature, they appear lower than those of similar objects studied 
so far such as MWC 480, HD 163296, and HD 34282 

- We found that change in UV opacity due to grain growth, 
and thus the resulting enhancement of the UV penetration, is 
a major process to explain the low CO column densities. 

- An important uncertainty in the PDR model resides in the 
estimate of the UV field intercepted by the disk. The amount 
of stellar UV photons is too small to produce the required CO 
photodissociation. The stars must be accreting to produce a 
sufficient UV excess. Since the accretion rates and the UV 
fields are both poorly known for most similar sources, this 
remains a limiting factor for chemical modeling of disks. 

- CQ Tau may be explained by a high UV flux and a normal 
g/d ratio. However, obtaining column densities just around 
10'^ cm"^ requires a rather fine balance between incident 
UV flux and UV penetration. 

- The lower gas temperature of MWC 758 is not well repro- 
duced by the PDR model. A reduced g/d ratio may be fa- 
vored for this source. Alternatively, the uncertainties in the 



overall thermal balance modeling may be the cause of the 
discrepancy. 

- Although the observed CO gas is too warm to allow sticking 
onto grains, we find that the temperature of large grains to 
be low enough to prevent CO from being released from the 
grain surfaces. The apparent low CO to dust ratio may thus 
be the result of the disk thermal history. 

- Finally, the inclination angle we measured in the outer disk 
of CQ Tau (29 ± 2") appears to be incompatible with the fact 
that this object is classified as an UX Ori star One possibil- 
ity left is that the inner disk is warped due to dynamical in- 
teractions with large bodies. Accurate determinations of the 
inclination of the inner and outer disks in more UXOr are 
required to progress in resolving this issue. 

Our study clearly indicates that a low apparent CO to dust ratio 
does not necessarily imply a low H2 to dust ratio, and thus that 
CO is not an unambiguous tracer of gas dispersal. 
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Appendix A: The Extinction Curve witli Grain 
Growth 

The extinction curves for several a+ can be seen on figure IaTI 
The peak around A = 2.2 lO^^m is due to small graphite grains 



(lMathisl[T994l) . the feature around A - I lO^^m is due to sil- 
icates. The curve for a+- 10"^ crosses the one for a+- 10"^ 
cm at 5.55 lO^^m, and the Ay extinction is calculated at A = 
5.5 lO^^m. This explains why the Nh/Av ratios for these two 
dust size distributions are so similar (table lATb- 



Table A.l. Extinction factor N///A,, 
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Nn/Ay ratio as a function of the maximum radius of the grain size dis- 
tribution for a gas to dust ratio of 100, a_= 3nm and y = -3.5. 
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Fig. A.l. Extinction curves as a function of the maximum radius 
of the grain size distribution for a gas to dust ratio of 100, a_= 
3nm and y - -3.5. 



When the grains are bigger the opacity in the mid-plane is 
smaller so that the photo-dissociation layer is more extended to- 
ward the mid-plane and CO is less abundant (Figure IA!2] |. In the 
case of relatively big grains («+= 10 fj.m) with a gas to dust ratio 
of 10 is the dominant form of carbon in the mid-plane for 
R > 200 AU. The radial distribution if the column density are 
presented on figure |A.3l 



Appendix B: Effects of UV field 



For smaller UV fields, CO is more abundant (figure IB.1IB.2| |. 
The most important part of the UV spectrum is the one be- 
tween 900 and 1200 A where the lines of H2 and CO are, be- 
cause the self-shielding of CO is important (see the effect for the 
two gas to dust ratio figure|4]l. We perform a calculation with a 
modified UV field to mimic the shape of the spectrum taken by 
[jonkJieid et all (|2007|) (a cut-off below 1200 A) with the standard 
scaling factor 10^. The H-H2 transition is dramatically affected: 
it occurs for a visual opacity less than lO--' magnitude instead of 
about 10-' for the x ^ ^0* and 10-^ magnitude with;t' =10-' for 
the Draine field shape. CO is less affected but the transition C- 
CO is also shifted towards the low opacity region. However the 
UV excess in Herbig Ae star is important and the spectrum shape 
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between 912 and 1600 A (where photodissociation takes place) 
is quite well represented by a scaled Draine field (see Figure|7]i. 

Appendix C: Turbulence 

We present here the effect of a multiplication and division by a 
factor 2 of the Doppler width of H2 lines on the surface density. 
The effect is maximum in the case of big grains with the ther- 
mal balance calculated but remain negligible for the C-bearing 
species (figure ICTT i. 
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Fig. A.2. Grain size effect on molecular abundances. Vertical distribution through the disk of the abundance of H, H2, C^, C and CO 
and of the temperature ( no thermal balance calculated) at the radii 100, 200 and 300 AU for the models with 0+= 0.1, 1, lOfim and 
1 mm. In the latter case the dust disk is so optically thin that we never reach A,, = 1 . 



Chapillon et al.,: Disks around CQ Tau and MWC 758 



17 



n 



thick : g/ d ratio = 1 
da£bed^l^g/d_ ratio = 100 



= 0-1 j^m 



100 



150 



200 
Radius (AU) 



250 



n 



thick : g/ d ratio = 1 
dashed : g/d ratio = 100 



= ^QfJ.m 



~r 



100 



150 



200 
Radius (AU) 



250 



y,= 1,e4 



CO 

C 



300 

. y,='',e'^. 



CO 

c 




300 





thick : g/d ratio = 10 
dq§he£j_; g/d ratio = 100 



100 



150 



200 
Radius (AU) 



250 



thick : g/d ratio = 10 
dc^hed_: _g/d ratio = 100 



100 



150 



200 
Radius (AU) 



250 



y.= 1.e4 




CO 

Z 



300 
. y.= 1.e4 



---tC 



CO 

c" 




300 
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Fig.B.l. UV field effect on abundances. Vertical distribution through the disk of the abundance of H, H2, C^, C and CO and of the 
temperature (no thermal balance calculated) at the radius 100, 200 and 300 AU for the models with small grains , a scale factor on 
the Draine field - 10^, 5.10^ and 10'* at 100 AU and with the modified UV shape (see text for explanation). 
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Fig.B.2. UV field effect on abundances. Vertical distribution through the disk of the abundance of H, H2, C^, C and CO and of the 
temperature (thermal balance calculated) at the radius 100, 200 and 300 AU for the models with big grains, a scale factor on the 
Draine field = lO"', 5 10^ and 10"* at 100 AU and with the modified UV shape (see text for explanation). Dust temperatures are 
plotted for the extreme grain sizes (a+and «_) and an intermediate value C^^" )- 
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Fig. B.3. Grain size effect on column densities for a weak UV field. Radial distribution of the column density of C^, C and CO for 
the model with weak UV field (x =10^ at 100 AU), a+- l/jm (left) and a+- 1mm (right). 




100 150 200 250 300 100 150 200 250 300 

Radius (AU) Radius (AU) 



Fig. B.4. Thermal balance effect on column densities. Radial distribution of the column density of C^, C and CO for the model with 
weak UV field (x = 10^ at 100 AU), a+= 1mm without (left) and with (right) the thermal balance calculated. 
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Fig. C.l. Turbulence effect over column density. Radial distribution of the column density of H, H2 (left), C^, C and CO (right) for 
three values of the Doppler width with the thermal balance calculated.a+= 1mm (bottom) 



